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ABSTRACT 

To investigate order-order interfaces, we perform multimagnetical Monte Carlo simu- 
lations of the ID and 3-D Ising model. Stringent tests of the numerical methods are 
performed by reproducing with high precision exact ID results. In the physically more 
interesting 3Z? case we estimate the amplitude of the critical interfacial tension. 
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1. Introduction 

Since long, there has been continuous interest in the properties of interfaces 
in Ising models. Past numerical studies were, however, hampered by a problem 
of principle. The surface tension per unit area F'^ between different states has a 
finite value. Thus in the canonical ensemble, where one samples with the Boltz- 
mann weights oc e~^^ , configurations containing interfaces with an area A are 
supressed by exponentiall factors e^^^ . Corresponding there is an exponentially 
fast increase of the tunneling time between pure phases when the system is simu- 
lated with localMontc Carlo (MC) algorithms. Recently this difficulty was overcome 
by the proposal to perform the MC simulation for a multimagnetical ensemble,^ a 
natural extension of the multicanonical ensemble introduced in (Ref. 2). The origi- 
nal canonical ensemble is then obtained by re-weighting,^'* while the slowing down 
becomes reduced to a power law close to ^ V^. 

For our models, spins = ±1 are defined on sites of a square lattice of volume 
V = with periodic boundary conditions and the symbol < i,j > is used to 
denote nearest neighbors. The partition function of the Ising model is given by 

Z = eM-PH). (1.1) 

configurations 

The Ising Hamiltonian H 

H = Hi - hM with i?/ = - ^ s^sj and M = ^ (1.2) 

<i,j> i 
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contains the nearest neighbor interaction term Hj and a term which couples the 
external magnetic field h to the magnetization. In our case is ft, = 0. In the broken 
region (3 > fic the magnetic probability densities Pl{M) are double peaked and 
Pl{M) = Pl{—M), as the model is globally Z{2) symmetric. The positions of the 
maxima are ztM^^^ and the distribution takes its minimum at p™'" = Pl(0). 

The surface tension is the free energy per unit area of the interface between 
coexisting phases. Here we consider order-order interfaces in the broken region 
f3 > Pc- The interface tension can be defined from the finite size scaling (FSS) 
of the ratio pmm/pmax^ which for large L takes the form^'^ 

pmin 

= const exp [-2L^-ip^(l + 0{L-'))] (1.3) 

on lattices with periodic boundary conditions. Physically this definition assumes 
that at values of the mean magnetization M = a rectangular domain, enclosed 
by two interfaces and spanning the lattice via the periodic boundary condition, is 
formed. The finite lattice interfacial tension is then defined as 

1 ^min 

n = ^^In;^, (1.4) 

and = lim^^^ Fl. So the final estimate of F^ requires a FSS extrapolation 
towards L = oo. Prom equation "(1.3)" one gets the fit: 

2. Multimagnetical Monte Carlo 

To get a representative sample of interfaces one has to generate many configu- 
rations with M = on large lattices. But, in a canonical simulation this is almost 
impossible. As an example see figure 1. The logarithmic scale displays that more 
than twenty orders of magnitude are involved, i. e. pmm^pmax ^ iQ-^o £qj. ^j^g 
largest lattice with L = 100. The standard MC algorithm would only sample con- 
figurations corresponding to -Pi^o" if one could generate of the order O(10^°) or more 
statistical independent configurations. 

Here we overcome this difficulty by sampling configurations with a multimag- 
netical weight factor 

P]^"(M) ^cxpiaL{M) + hLiM).dM - l3Hi) for -V < M <V (2.1) 

instead of sampling with the Boltzmann factor = exp(— /3if/). The function 
hL{M) defines an M-dependent effective external magnetic field such that the re- 
sulting multimagnetical probability density is flat for — M™^^ < M < M™"*^. We 
chose the ansatz: 

, , _ J for M < -Mf and M > M™'^'^; 
hL[M) - |2-i/3-MlogPL(M)-logPL(M + 2)] elsewhere. ^ > 



Fig. 1. Boltzmanii probability distribution Pi^[M) of the inagnctizatioii for different lattice sizes. 



The function ai,(M) is then defined by the recursion relation 

aL{M + 2) = aL{M)+p[hL{M)-hLiM + 2)] (M + 2), ^^(-M^''") = 0. (2.3) 

Once hL{M) is given, ai(M) follows automatically. With this choice the resulting 
multimagnetical probability density will be approximately flat: 

P™"(M) = n{M) ■P™'"(M) w const. (2.4) 

Here n{M) is the magnetical density of of states (for fixed temperature The 
standard Markov process is well-suited to generate configurations which are in equi- 
librium with respect to this multimagnetical distribution. The canonical probability 
density Pl{M) is obtained from P]^"^{M) by re-weighting^'^: 

Pl{M) = c PE'"'{M)eM-aL{M) - hL{M)pM). (2.5) 

The constant c is obtained by imposing the appropriate normalization on Pl(M). 

For small systems Pl{M) can be calculated by performing standard MC sim- 
ulations, and hL{M) follows directly from "Eq. 2.2". On larger systems wc get 
hL{M) by making every time a FSS prediction of Pl{M) from the already con- 
trolled smaller systems. 

To compare the efficiency of our method with standard MC we measured for 
the 2D Ising model with (3 = 0.5 the tunneling time t£. As before^'^ we define 
the tunneling time as the average number of updates needed to get from a 



Fig. 2. Tunneling times versus lattice size. The upper data are for the standard MC and 
the broken line extrapolates by means of a fit into the region where no data exist. The lower data 
points are obtained with our multimagnetical MC. 

configuration with M = M'"'" = -M'^'''^ to a configuration with M = M^'^'' 
and back. In figure 2 we display on a log-log scale both the tunneling times for the 
multimagnetical MC and the heatbath algorithm. While there is an exponentiall fast 
increase of for the heatbath algorithm, the increase of is for the multimagnetic 
MC of the type of a power law divergence. The ratio 



is a direct measure for the improvement due to our method. R increases from a 
factor 4 for the smallest lattice {L = 2) up to i? ~ 450 for L = 16, the largest lattice 
size where it was with our statistics possible to get data from standard MC. The 
extrapolation to L = 100 yields ii w 6.1 x 10^^, i.e. an improvement by more than 
fifteen orders of magnitude. 

3. Numerical Results 

For the two dimensional Ising model we performed multimagnetical simulations 
at the critical temperature (3 = l3c = ln(l + \/2)/2 = 0.44068..., at /3 = 0.47 and 
0.5 with at least 4 x 10^ sweeps per lattice size. In each run additional 200, 000 
initial sweeps without measurements were performed for reaching equilibrium with 
respect to the multimagnetical distribution. We compared our L = oo estimates of 
the tension with the exact values which follow from Onsager's equation^ 



R — (heatbath) /r|,(multicanonical) 



(2.6) 



F' = 



2/3 -In 



+ 

1 _ e-2/3 



(/3 > /3c). 



(3.1) 



Fig. 3. ID eflfective tensions. The lines are our fits 



Both are collected in table 1. Good agreement is found in all cases. 

Table 1. Surface tensions for the ID Ising model 



/3 


mm 


exact 


/3c 


0.00033(34) 





0.47 


0.11526(80) 


0.11492... 


0.5 


0.2281(5) 


0.22806... 



In the case of the 3D Ising model we performed simulations at /? = 0.227, 
0.232 and 0.2439 with at least 4 x 10^ sweeps for every lattice size. As in the 
2D case 200,000 additional, initial swoops were performed in oacli run for roaohing 
equilibrium. In figure 4 wc display the effective tensions as functions of the lattice 
size together with asymptotic fits. We notice that finite size effects play a more 
important role in three dimensions than two and are more complicated, too. The 
non-monotone behaviour shows that it is necessary to use large enough lattices to 
estimate the interface tension. 

The value /3 = 0.232 was chosen, because it enables an comparison with the 
recent literature. In (Ref. 8) cluster improved estimators were used to calculate 
correlations for L = 8 — 14 in a cylindrical geometry. Fitting the obtained tunneling 
mass gaps yields F'^ = 0.03034 ± 0.00015 for the surface tension, when systematic 
errors are admitted, this is consistent with our value. All our F^ (see table 2) values 



Fig. 4. 3D effective tensions. The lines are our fits. 



are much higher that the old estimates of (Ref. 6) , which had to rely on far too small 
lattices. Our estimates are also consistent with the valueses presented by Miinster^ 
and Potvin^° in other talks on this workshop. 

Table 2. Surface tensions and the amplitude Fq for the 3D Ising model 



P F^ F^ 



0.227 0.01293(17) 1.455(9) 
0.232 0.03140(14) 1.581(7) 
0.2439 0.07403(30) 1.513(6) 



Of particular interest in the 3D Ising model are amplitude ratios^^^^^ which in- 
volve the amplitude Fq of the critical 3D interfacial tension, because they can be 
compared with a variety of experimental results^'^^'^^ for fluids which are supposed 
to populate the Ising model universality class. The other amplitudes are fairly 
accurately known, and the uncertainty in the ratios stems mainly from Fq. Old 
estimates of Fq (Ref. 6) led to amplitude ratios which are in disagreement with 
the experimental results"^^'"*^^, whereas Mon obtained consistency by calculating the 
excess energy between a system with an interface imposed by an antiperiodic bound- 
ary condition and a system without such an interface and with periodic boundary 
conditions. 



The asymptotic behaviour of the interface tension for (3 ^ l3c is 



F' = F^t", (3.2) 

where t = {1 — Pc/P) is the reduced temperature and /x = 1.26 in 3D. Our t-values 
are small enough to apply this formula. The results for the amplitude are included 
in table 2. By averaging over our three values we get 

F^ = 1.52 ±0.05, (3.3) 

compared to Mon's value F^ = 1.58 ± 0.05 (Ref. 15). 

4. Conclusion 

Multimagnetical simulations allow to study the magnetic probability density in 
the broken region with a hitherto unreached precision. Our numerical calculations 
for the 2D Ising model agree well with the exact results. For the 3D Ising model 
we obtain new surface tension estimates,which agree resonable with (Ref. 8,9,10), 
and for our amplitude good agreement is found with Mon (Ref. 15). 
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